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O ■ Abstract 

Evidence from the WMAP polarization data indicates that the Universe may have been 
reionized at very high redshift. It is often suggested that the ionizing UV flux originates 
from an early population of massive or very massive stars. Depending on their mass, such 
stars can explode either as type II supernovae or pair-instability supernovae, or may entirely 
collapse into a black hole. The resulting neutrino emission can be quite different in each case. 
We consider here the relic neutrino background produced by an early burst of Population III 

a ■ 

stars coupled with a normal mode of star formation at lower redshift. The computation is 
performed in the framework of hierarchical structure formation and is based on cosmic star 
formation histories constrained to reproduce the observed star formation rate at redshift 
z < 6, the observed chemical abundances in damped Lyman alpha absorbers and in the 
intergalactic medium, and to allow for an early reionization of the Universe at z ~ 10—20. We 
find that although the high redshift burst of Population III stars does lead to an appreciable 
flux of neutrinos at relatively low energy [E v rs 1 MeV), the observable neutrino flux is 
dominated by the normal mode of star formation. We also find that predicted fluxes are 
at the present level of the SuperK limit. As a consequence, the supernova relic neutrino 
background has a direct impact on models of chemical evolution and/or supernova dynamics. 



1 Introduction 



One of the most surprising results contained in the first year data obtained by WMAP PQ 
was large optical depth implying that the Universe became reionized at high redshift in the 
range, 11 < z < 30 at 95% CL. This is much higher than previously believed, and it has 
been proposed that a generation of very massive stars is necessarily responsible 0. Even 
a brief period of massive star formation at high redshift would have notable consequences 
for chemical evolution, in particular, the metal enrichment of the interstellar medium (ISM) 
and intergalactic medium (IGM) jHllUEl- When combined with the observed cosmic star 
formation rate (SFR) at z < 6 [6 , a coherent picture of the star formation history of 
the Universe unfolds: the first stars (Population III) described by a top-heavy initial mass 
function (IMF) were formed in primordial, metal-free structures with masses of order 10 7 M Q . 
Once the metallicity achieved a certain critical level (of order 10~ 4 times the solar metallicity 
[Zj), the massive mode of star formation yielded to a more normal distribution of stellar 
masses at a rate over an order of magnitude larger than the current star formation rate. The 
cosmic SFR is observed to have peaked at redshift z ~ 3. 

Among the consequences of this newly emerging view of star formation is the predicted 
enhancement in the rate of core collapse supernovae. In addition to the sharp spike of 
supernovae at very high redshift due to the explosions of stars responsible for the early 
epoch of reionization, the enhanced SFR of the normal mode of star formation at redshifts 
z < 6 leads to a supernova rate which is approximately a factor of 30 times the current rate, 
and a factor of 5 times the observed rate at z ~ 0.7 [HJ. 

Another consequence of an enhanced SFR and SN rate is the resultant neutrino back- 
ground spectrum produced by the accumulated core collapse supernovae. Early estimates 
of the (anti-e)-neutrino flux based on simple models of Galactic chemical evolution were in 
the range of 1 - 10 cm~ 2 s~ 1 MeV _1 pj, as were later estimates ^Hl El E] based on peak 
SFRs of order 0.1 M yr _1 Mpc _3 corresponding to supernova rates of order 10 -3 yr _1 Mpc~ 3 . 
In PHI, the first attempt at obtaining the neutrino flux from very massive Population III 
stars found fluxes at levels exceeding 10 cm _2 s _1 MeV _1 , though the peak occurs at lower 
energy due to the redshifted spectrum. Current estimates of the SFR peak at values 3 times 
higher at a redshift z ~ 3 [6 j with correspondingly higher SN rates. More startling is the 
possibility that SN rates at high redshift (z ~ 15) may be as large as 6 x 10 -2 yr~ 1 Mpc -3 . 
Here, we incorporate fully developed chemical evolution models which trace the history of 
pre-galactic structures as well as the IGM and are based on a ACDM cosmology and include 
a Press-Schecter model [13] of hierarchical structure formation. We adopt the approach 
and chemical evolution models of Daigne et al. 01 |5] and consider several bimodal IMFs, 
each with a normal component of star formation as well as a massive component describing 
Population III stars. Given an IMF and a respective SFR, one can calculate the chemical 
history of cosmic structures and the IGM, the reionization efficiency, and as we do here, the 
expected supernova relic neutrino (SRN) flux. 

In section |2] we describe our method for calculating the SRN flux, and in section |S] we 
discuss the star formation models considered. In section 0] we present our results for the SRN 
flux in a variety of models describing Population III stars. Opportunities for SRN detection 
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are examined in section and the effects of neutrino oscillations are discussed in section El 
Our conclusions are given in section [7| 



2 Supernova Relic Neutrino Background 

In all of the models we will consider, star formation begins at high redshift, dominated 
initially by massive stars which may explode as core collapse or pair-instability supernovae 
and provide for the reionization of the IGM. Each explosion, regardless of type, releases 
most of the star's gravitational energy in the form of neutrinos with a specified energy 
spectrum and flux. Given a chemical evolutionary model, or more specifically, given a rate 
of star formation along with the IMF, the integrated contribution of SN to the neutrino 
background can be computed. 

The expected differential flux of neutrinos at Earth with energy E can be expressed as 
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where N l/a (m) is the total number of neutrinos of a given species, a, emitted in the core 
collapse of a star of mass m, is the normalized IMF, ip is the SFR per unit comoving 
volume, r(m) is the lifetime of a star of mass m, and ^7 represents the neutrino spectra in 
the comoving volume at energy E' = E(l + z). The integration limits M m j n and M max are 
the minimum and maximum masses in each model for which supernovae occur, and Z{ is the 
initial redshift over which star formation occurs. It is assumed that star formation continues 
to the present [z — 0). 

When a star undergoes core collapse, the mass of the remnant is determined by the mass 
of the progenitor. We assume that all stars of mass m > 8 M will die as supernovae. For 
stars of mass 8 M Q < m < 30 M , the remnant after core collapse will be a neutron star 
of m Pd 1.5 M Q . More massive stars fall into two categories; black holes and pair instability 
supernovae. Pair instability supernovae are thought to occur for stars with 140 M < m < 
260 M Q , in which case the explosion leaves no remnant. All other stars collapse to form 
black holes. Stars with 30 M < m < 100 M become black holes with mass approximately 
that of the star's helium core before collapse [15J. We take the mass of the Helium core to 
be 

13 

M He = --(m-20M e ) (2) 

for a star with main sequence mass m [16 a . We assume that stars with m > 260 M collapse 
entirely to black holes. 

The energy emitted in each core collapse, E cc corresponds to the change in gravitational 
energy, 99% of which is emitted as neutrinos J7j. In the cases where collapse results in a 
neutron star, E cc = 5 x 10 53 ergs. For stars that collapse to black holes, E cc is proportional 
to the mass of the black hole. Pair instability supernovae experience a much more powerful 
explosion than core collapse supernovae, however few neutrinos are emitted and with very low 
energies such that they would not be observed [IB] . We address this issue further in sectional 



2 



Although several studies find a distinct hierarchy in the partitioning of neutrino luminosity 
among the species during the different luminosity phases of core collapse, equipartition of 
the total energy emitted by the star is generally accepted El C2| For a comparison of 
luminosity hierarchies found in recent simulations, see Keil et al. [TH] . 

In black hole formation, the neutrino luminosity is nearly constant for the first few sec- 
onds until the event horizon overtakes the neutrinospheres. Once the neutrinospheres are 
inside the event horizon, the luminosity consists of neutrinos with lower average energies 
escaping from the outer layers of the star [201 E2] • We assume that each electron neutrino 
carries an average energy (E Ue ) = 13.3 MeV, which is a reasonable approximation for the 
average neutrino energy over the two luminosity phases of core collapse to a black hole. 
For supernovae which do not collapse to black holes, this energy is consistent with recent 
simulation data ^H]- The charged current reactions that prevent neutrinos from emerging 
from the star are z/en — > pe~ and z/ep — > ne + . The different trapping reactions result in 
different neutrinosphere radii, and therefore different average energies for v e and v e . We as- 
sume (E Pe ) = 15.3 MeV which is the average energy over the two luminosity phases as above, 
following The other species, denoted u x , undergo only neutral current interactions. The 
mechanism that governs their average temperature at emission is more complicated, but the 
generally accepted hierarchy is (E Ue ) < (E Pe ) < (E Vx ). We have taken (E Vx ) = 20 MeV. The 
total number of v a emitted by a star during core collapse is given by 

N v =-^t- (3) 

We discuss the sensitivity of our results to these choices in section El below. 

The neutrino spectra can be described by a normalized Fermi-Dirac distribution, 

dP a 2 E' 2 

dE> 3C 3 T3 e E '/ T « + 1 U 

where T a = 180( 3 (E Ua ) /77r 4 is the effective neutrino temperature taken to be independent 
of the mass of the star. We assume a flat ACDM cosmology with 
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where fi A = 0.73, Vt m = 0.27, and h = 0.71 pp. 



3 Star Formation Models 

The cosmic star formation histories we consider here have been adopted from detailed chem- 
ical evolution models [3J E|. These models are bimodal and are described by a birthrate 
function of the form 

B(m,t) = 0i(m)V>i(*) + <h( m )M t ) ( 6 ) 

where <f>i(2) is the IMF of the normal (massive) component of star formation, and 1/^1(2) 
is the respective star formation rate. The normal component contains stars with masses 
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between 0.1 and 100 M and is primarily constrained by observations at low redshift (z < 6). 
The massive component operates at high redshift and is required by the evidence for the 
reionization of the Universe at z ~ 17. Both components can contribute to the chemical 
enrichment of galaxy forming structures and the IGM. We consider three different models 
of the massive mode as described below. 

Given an IMF and SFR, it is straightforward to compute the rate of core collapse super- 
novae, 

SNR= / dm <t>(m)il> (t - r(m)) , (7) 

J max(8MQ,m ul i n (t)) 

where m m m (i) is the minimum mass with lifetime less than t. Then for each model, the 
SRN flux is calculated using eq. ((IJ). In the two following subsections, we describe the two 
evolution models we have used, and their corresponding SFR and IMF. 



3.1 Stage 1: Simplified Models 

We begin by describing a simplified set of hierarchical models jlj. We assume star formation 
to occur between Zi = 20 and Zf = 0. The IMF for each mode is similar to a Salpeter mass 
function, 

0i(m) oc m- (1+x) (8) 
but has a slightly steeper slope of x = 1.7. Each IMF is normalized independently by 

sup 

dm mcj)i(m) = 1 , (9) 

differing only in the specific mass range of each model. 

The mass in cosmic structures will be denoted M struct (t), which includes both the mass 
in stars and the mass in gas of the ISM. The mass in the IGM is Mjqm and the total mass 
M tot = Mjqm + M struct is of course constant. The mass in structures evolves as 

= a b (t) - o{t) , (10) 

where a&(i) is the cosmic baryon accretion rate due to the structure formation process and 
o(t) is the outflow from structures. In our simplified set of models, which we denote as stage 
1, we will ignore the effects of outflow and set o(t) = 0. We will also assume structure 
formation to be exponentially decreasing from time t — 0, 

a b (t) = -Mtote-^' , (11) 

where a = 0.1 is the fraction of the total mass which is eventually accreted by structures 
and Ts is the timescale of the accretion process. Here we consider both r s = 0.01 Gyr and 
t s = 0.2 Gyr We also assume that star formation begins when the baryon fraction in 
structures is 1%. This is our initial condition at z = 20 and corresponds to an estimate 
of the minimum baryon fraction where sufficient dissipation occurs to allow star formation 
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In our more complex models we retain the condition on the onset of star formation, but 
determine its redshift using a detailed model of hierarchical structure formation (see next 
subsection) . 

The evolution of the gas mass in the structures is given by : 

= + e (t) + a h (t) - o(t) . (12) 

at 

The stellar mass in the structures is simply M strU ct — Mjsm- Each term in equation (JT2|) 
accounts for part of the gas budget of the ISM. The first term corresponds to the loss of gas 
through star formation while the second term corresponds to the ejected gas when the star 
dies. Here, we will use the instantaneous recycling approximation (IRA) to evaluate the rate 
at which gas is returned to the ISM, e(t) 



e(t) = dm 4>(m)ip(t)(m - m r ) = Rt/j(t) . (13) 

io.9M 

where m r is the mass of the remnant, and R is the IRA return mass fraction. 

The normal mode of star formation in [3] is referred to as Model and provides a standard 
star formation history, with stellar masses in the range 0.1 M < m < 100 M . The SFR 
for Model is proportional to the gas mass fraction in cosmic structures, a = MisM/M struct , 

in = v x a(t) , (14) 

where v x = M struct (t) /r\ with T\ = 5 Gyr which is a typical timescale for star formation in 
the galactic disk. This model alone is inadequate for high redshift reionization. 

We consider three different models, labeled Models 1, 2a, and 2b to describe the massive 
mode. They are distinguished by their respective stellar mass ranges. In Model 1, the IMF 
is defined for stars with masses, 40 M < m < 100 M . All of these stars die in core collapse 
supernovae leaving a black hole remnant. They all contribute to the chemical enrichment of 
the ISM (and IGM when o(t) ^ 0). This period of star formation is brief and is described 
by a SFR of the form 

^ 2 = v 2 e- t/T2 (15) 

where vi = f 2 M strU ct{t) /t~2 with a characteristic timescale r 2 = 50Myr. The constant 
/ 2 = 4.5%. 

Model 2a is described by very massive stars which become pair instability supernovae. 
The IMF is defined for 140 M < m < 260 M and the SFR for this model is the same as that 
for Model 1, but must be reduced by a factor of 8 due to constraints on metal abundances 
in the ISM. We consider the best case scenario for observation, where the energy emitted in 
neutrinos is the same as that for ordinary core collapse supernovae, but the average neutrino 
energy is (E Pe ) = 1.2 MeV HE]. 

The most massive stars are considered in Model 2b and fall in the range 270 M < m < 
500 M , with the SFR as in Model 1. These stars entirely collapse into black holes and 
do not contribute to the chemical enrichment of either the ISM or IGM. Collectively we 
will refer to stage 1 models as 1.0, 1.1, 1.2a, and 1.2b. Unless otherwise noted, Models 1.1, 
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1.2a, and 1.2b (as well 2.1, 2.2a, and 2.2b denned below) will correspond to the full bimodal 
model, i.e., they include the normal mode of star formation as well as the particular massive 
mode. 



3.2 Stage 2: Hierarchical Models 

We also consider a set of models with a more sophisticated treatment of the hierarchical 
growth of structure which we will call here stage 2. Complete details of this model can be 
found in [Sj. Rather than assuming an analytic form for the formation of structure as in eq. 
(jllj) . we use the Press-Schechter formalism and take 
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where fb is the baryon fraction in structures 



fb, struct ('') 



J Q °°dM Mf PS {M,z) 



(17) 



In eq. (JHj), fps(M, z) is the distribution function of halos and is computed using the method 
described in , and M min is the minimum halo mass in which stars are formed. As before, 
we require an initial baryon fraction fb = 1%, which for example for M mui = 10 7 M 
corresponds to an initial redshift for star formation, z = 16. 

In jH], it was found that better fit to the global star formation rate and supernova rates 
is obtained when a normal mode SFR of an exponential form 
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(18) 



is used, which corresponds to a SFR dominated by elliptical galaxies. Best fits for vi and T\ 
are given in Table 1. In these models the outflow is non-zero and the details for computing 
the outflow are given in |EJ E] . The overall efficiency of outflow is parameterized by e whose 
value is also given in Table 1. The instantaneous recycling approximation is no longer used 
here, and ejection rates depend on stellar lifetimes, r(m). This amounts to replacing ijj(t) 
with ip(t — r{m)) in the integral in eq. (|13|) . 

The SFRs for the massive modes in stage 2 are determined by the metallicity in the ISM, 



V>2 



-z/z c 



(19) 



where Z cr - lt = 10~ 4 Z is the critical metallicity at which Population III star formation ends 
[7j. The IMF of both modes in this case has a slope £i( 2 ) = 1.3. The massive mode SFR 
parameters are given in Table 2 for all models considered. The mass ranges and neutrino 
average energies in stage 2 are the same as those in stage 1 respectively. 
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Normal mode 





2qnit 


e 




T\ 


(M ) 






(Gyr v ) 


(Gyr) 


10 6 


18.2 


2 x 10~ 3 


0.2 


2.8 


10 7 


16.0 


3 x 10~ 3 


0.2 


2.8 


10 8 


13.7 


5 x 10~ 3 


0.2 


2.8 


10 9 


11.3 


io- 2 


0.2 


3.0 


10 11 


6.57 


1.5 x IO" 2 


0.5 


2.2 



Table 1: The model parameters for the normal mode of star formation (Model 2.0). Column 
1 indicates the input value of the minimum mass for star forming structures. Column 2 
is derived from column 1, having assumed that /& = 1% when star formation begins. In 
columns 3, 4, and 5, parameter values for the efficiency of outflow and the SFR are given. 
The slope of the IMF is x = 1.3 for all models. 



Massive mode 



Model 


M min 


e 




V2 




(Mo) 






(Gyr" 1 ) 


2.1 


10 7 


2 x 10" 


-3 


60 


2.1e 


10 7 


6 x 10" 


-5 


340 


2.2a 


10 7 


1 x 10" 


-3 


9 


2.2ae 


10 7 


8 x 10" 


-5 


40 


2.2b 


10 v 


3 x 10" 


-3 


100 



Table 2: Parameter values for the massive starburst Models 2.1, 2.2a and 2.2b. Column 1 
indicates the model number and column 2 the input value of the minimum mass for star 
forming structures. In columns 3 and 4, we show the adopted outflow efficiency and massive 
mode SFR. 

When a massive mode is added to the normal mode described by Model 2.0, the outflow 
efficiency must be adjusted so as to avoid the overproduction of metals in the IGM. However, 
there is a degeneracy in the massive mode parameters e and z/ 2 . In models labeled 2.1 and 
2.2a, the massive mode contributes roughly 50 % of the IGM metallicity at a redshift z = 2.5. 
By increasing vi and decreasing e, this contribution can be increased to 90 % and at the 
same time increases the ionization capacity of the model. These cases are labeled 2.1e and 
2.2ae 1 . 

1 Since stars associated with Model 2b do not contribute to element enrichment, there is no Model 2be. 
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4 Results 



4.1 Neutrino Background 

Stage 1 and Stage 2 models are discussed individually in sections 14.1.11 and 14.1.21 Then 
comparisons are made in section 14. 1 .HI with other results. We will for the most part consider 
only u e , as they are the most easily detected at water Cerenkov detectors. We will return to 
v e when we address thermonuclear neutrinos and detection. 

4.1.1 Stage 1 

We begin by examining the neutrino production from simplified models, labeled here as 
stage 1. As described above, these models employ the IRA and as such can to a large extent 
be treated (semi)-analytically. In this class of models, outflows are ignored and hierarchical 
growth is simply modeled by eqs. (jlOj) and (jlljh All of the Population III models we consider 
are bimodal (see eq. (JBJ)) and combine a normal stellar distribution (model 0) along with 
a massive mode. As a reference point, we first compute the expected neutrino flux from 
the normal mode alone. This is shown in Fig. El where we show the flux of v e only. The 
peak flux is about 2 cm" - 2 s _1 MeV _1 and occurs at E v ~ 2 MeV. While this flux is small 
compared to solar neutrino fluxes (for example, the flux of 8 B neutrinos from the sun is of 
order 10 6 cm _2 s _1 MeV~ 1 ), it is large compared to the non-localized atmospheric neutrino 
flux which is less than 10 -2 cm _2 s _1 MeV _1 (see e.g. [TU]). 

Fig. ^ shows the resultant neutrino flux for two choices of the baryon fraction in struc- 
tures, M struct I 'M to t = -1 (large dashes) and for comparison, M struct /M tot = .01 (small dashes) 
when the growth of structures is neglected (i.e. a,b = 0). As expected, the neutrino flux 
is in direct proportion to the baryon fraction. We also show in Fig. ^ the neutrino flux 
when a;, 7^ 0, and the mass of the structure grows. We show results for two different growth 
constants, t s = 0.01 and 0.2 Gyr. Since the final baryon fraction is 10% in each case, we see 
that the integrated neutrino flux is quite similar. 

The individual fluxes for the massive modes of Models 1.1, 1.2a, and 1.2b are displayed 
in Figure El For each of the massive modes, the duration of star formation is very brief, as 
the SFR is characterized by a time constant T2 = 0.05 Gyr. The neutrino fluxes for Models 
1.1 and 1.2b, shown by solid (blue) and dotted (black) curves respectively, are significantly 
larger than those found for a normal population of stars. In these cases, the fluxes are 
approximately 30 and 95 cm _2 s _1 MeV _1 , though the peak of the spectrum occurs at lower 
energy, E v ~ 0.6 MeV. We note that the SRN flux from Model 1.2b is very similar to the 
flux due to the population of rotating 300 M stars with a SFR peaked at z = 17 in Iocco 
et al. PS]- The peak height obtained in our calculations is larger due primarily to the fact 
that the integrated fraction of baryonic matter in population III stars is about four times 
greater in our model. It has been noted that if the collapse to a black hole proceeds without 
rotation, the neutrino luminosity will be diminished by ~ 2 orders of magnitude |20j . so this 
is really an upper limit to the flux assuming rotation. 

The flux from Model 1.2a is much smaller than the fluxes from Models 1.1 and 1.2b, and 
peaks at a lower energy. This is a consequence of less energy being released in neutrinos by 
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Figure 1: Electron anti- neutrino fluxes from Model 1.0 with different baryon accretion rates. 
We have chosen two models with a& = with different baryon fractions, M st ruct/M t ot = -1 
(red dashed) and M struct /M tot = .01 (green dotted). We have also chosen two models with 
structure growth using r s = O.OlGyr (thick black), and r s = 0.2 Gyr (solid blue). In the 
latter cases, a = 0.1. 

pair instability supernovae than by supernovae that collapse to form black holes, and that 
the average energy of each neutrino is limited by silicon burning and photodisintegration |18j . 

Comparison of Figs. ^ and EJ shows that the normal mode, though smaller in its peak 
flux, dominates the full bimodal spectrum at energies E > 2.5 MeV. Recall that the massive 
mode is very localized at high redshift. As a result, neutrinos produced by the massive mode 
have energies which are redshifted from their initial value. In contrast, the normal mode of 
star formation is peaked at z ~ 2 — 3, which produces a broader spectrum today. 

Finally, Figure H3 shows the total flux for the bimodal IMFs in the case that the massive 
component is either Model 1.1, 1.2a, or 1.2b. At low energies, the flux is dominated by the 
massive component. The insert shows the low energy peak due to PISN. 

4.1.2 Stage 2 

We now use the full numerical results of the models described in jS] . The overall form of these 
models is similar to those in stage 1, i.e. they are described by a bimodal birthrate function 
of the form in eq. ©. The mass ranges for models 0, 1, 2a, and 2b are unchanged, though 
the slope of the IMF in each case is now 1.3. As described above, these models include both 
outflow and the hierarchical growth of structure and we no longer employ the IRA in any of 
our calculations. The SFR is shown in Fig. 0]for Models 2.1, 2.2a, and 2.2b as indicated by 
their respective mass ranges. In Models 2.1 and 2.2a, the massive mode contributes roughly 
equally with the normal mode to the IGM metallicity at z ~ 3. Since Model 2.2b produces 
no heavy elements, the SFR is chosen to be maximal to enhance its ionization potential. In 
each case, star formation begins at z ~ 16 when the baryon fraction in a structure of total 
mass 10 7 M is ft, = 1%. The massive burst ends when the ISM metallicity has reached a 



9 



0.5 1 1.5 2 2.5 3 
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Figure 2: Fluxes from the massive modes of Models 1.1 (solid-blue), 1.2a (dashed-red), and 
1.2b (dotted-black). 

critical value taken to be 10 -4 Z . The duration of the burst in Model 2.2b is somewhat 
more prolonged as the metallicity is produced solely by the normal component. 

Model 2.0 fluxes are plotted in FigureEl Results are shown for several choices of minimum 
halo masses, M min . The neutrino fluxes found here are roughly a factor of 10 times larger 
than that found in the semi-analytical model discussed in the previous subsection. This is 
partly due to the flatter IMF chosen here which greatly enhances the numbers of massive 
stars and hence the supernova rate and neutrino flux. As we will see below these fluxes are 
large enough to be probed by current detectors. As M min is increased, star formation occurs 
at later redshift and as a result, the peak of the neutrino flux is shifted slightly to higher 
energy. 

The massive modes of Model 2.1 and 2.1e fluxes are plotted in Figure |H1 for the specific 
choice of M min = 10 7 M Q which is the preferred case in jSJ. As seen in Figure HJ massive stars 
associated with Population III turn on at a redshift of approximately 16, but the duration of 
the burst is relatively brief. As a result, the peak of the flux distribution is at relatively low 
energy. More importantly, because of the brevity of the burst, the entire neutrino spectrum 
is redshifted down, in contrast to the Model 2.0 spectrum which extends to higher energy 
due to stars produced at lower redshifts. As expected, the more extreme model, 2.1e, has a 
peak flux which is about 5 times that found for Model 2.1. This is directly related to the 
increased SFR in Model 2.1e as characterized by the increase in v 2 - 

Similarly, we show in Figure the resulting flux from very massive Population III stars 
corresponding to Models 2.2a, 2.2ae and 2.2b. As before, the fluxes from Models 2.2a 
are relatively small and peak at very low energy as seen in the insert to the figure. In 
Figures |H1 © ) we show the total fluxes in Models 2.1, 2.2a, and 2.2b (2.1e and 2.2ae) with 
M min = 10 7 M . As one expects, the low energy spectrum is dominated by neutrinos 
produced in the massive mode, whereas the spectrum at higher energies (E u > 3 MeV), is 
indistinguishable between the models and dominated by the normal mode. 
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Figure 3: Total fluxes from Models 1.1 (solid-blue), 1.2a (dashed-red), and 1.2b (dotted- 
black). In all three cases, a?, ^ and r s = 0.01 Gyr. The insert shows the low energy peak 
due to PISN. 

4.1.3 Comparisons to other models 

The differences in the fluxes between stages 1 and 2 models can be attributed entirely to 
differing IMF slopes, the baryon fraction, and treatment of the SFR. In order to compare 
our calculations to previous ones, it is necessary to single out individual models and discuss 
the differences. 

Our Model can be compared to models described in jHJ EH HH E2] • In M the calculated 
neutrino flux peaks at ~ 8 cm _2 s _1 MeV _1 with neutrino energies of ~ 3 MeV. In the models 
considered, star formation only occurs for z < 5 and the average neutrino energy at emission 
was assumed to be a step function in which energies are generally a few MeV larger than 
our adopted value of (E Pe ) = 15.3 MeV. Since the flux changes as 



where x depends on z, this both reduces the peak height and shifts the peak to larger neutrino 
energy at detection. In addition, they include only supernovae that collapse to neutron stars, 
with 8 M Q < m < 50 M Q . In the absence of any component of stars which collapse to black 
holes, the flux should be smaller than that obtained by our calculations. Ando, Sato and 
Totani use a crude supernova model, but provide a more careful treatment of oscillations. 
As will be shown in section our results are qualitatively similar despite the fact that our 
integrated flux differs from theirs by an order of magnitude. 

The more modern approach taken in ^01 EI] is based on the cosmic star formation rate. 
However, our Model 2.0 spectra have significantly higher peaks than those in ^UJ [n]. First, 
improved determinations of the SFR at high redshift are higher by a factor of about 2. A 
more significant difference between our results can be traced to our inclusion of stars that 
collapse to black holes. The change in gravitational energy is much greater when the remnant 
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Figure 4: The SFR for Model 2.1 (solid-blue), 2.2a (dashed-red), and 2.2b (dotted-black) 
as a function of redshift. Each model is indicated by the mass range associated with the 
massive mode. In all three cases, M m ; n = 10 M . Data are taken from [Bj. 



is a black hole than when it is a neutron star leading to a larger value for E cc , and therefore 
the neutrino luminosity is also much greater. For comparison, if we consider only stars that 
collapse to neutron stars, our flux is reduced by a factor of ~ 5. 

In [12], the SFR is parametrized as a broken power law, flat for z > 1 and fit to the 
observed cosmic star formation rate. As in [HIED!, star formation begins at z = 5. The effect 
of neglecting star formation at higher redshift is inconsequential for the higher energy tail 
of the flux, but both diminishes the lower energy flux and shifts the peak forward. Both of 
these effects are consequences of the fact that neutrinos from the earliest supernova events 
will arrive at Earth with the lowest energies as a result of redshift. 

The flux from a massive mode of star formation was computed in [13 . Population III 
stars were assumed to be rotating 300 M stars with a SFR approximated by a delta function 
peaked at z — 17. We observe the same peak location and spectral shape in our Models 
1 and 2b. Our fluxes peak at larger values than theirs for two main reasons. First, our 
integrated fraction of baryonic matter in Pop III stars is larger. They use a fixed value of 
10~ 3 , whereas our value is derived for each model; 2.2 x 10~ 3 for Model 2.1 and 7.3 x 10~ 3 
for Model 2.2b. Second, in Model 2.2b the stars collapse entirely to black holes, while in 
their model the energy emitted in neutrinos corresponds roughly to the gravitational energy 
released when the stars collapse to black holes with masses equivalent to the mass of the 
helium core before collapse. As a result, a 300 M Q star in our model emits ~ 2.7 times the 
number of neutrinos as the same star in their model. 

4.2 Thermonuclear neutrinos 

In this section we consider briefly the neutrinos emitted by stars during the hydrogen-burning 
phase. The CNO cycle is the dominant neutrino production mechanism, responsible for 90% 
of the thermonuclear neutrino flux. The remaining 10% is due to the pp chains. To simplify 
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Figure 5: Fluxes from Model 2.0 for five choices of M min = 10 6 (solid blue), 10 7 (thick black), 
10 8 (dotted red), 10 9 (dashed green) and 10 11 (dot-dashed grey) M . 
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Figure 6: Fluxes from the massive modes of Model 2.1 (blue) and 2.1e (cyan). 



the calculations, we consider only the reactions given in Table below. 

The total number of neutrinos emitted by a star of mass m during hydrogen burning has 
been estimated to be [T3] 

Tfl 

N therm— nucl ~ 0.2 , (21) 

m N 

where is the mass of a nucleon, and the energy spectrum at emission is taken from 
Ref. |23 • The total flux of thermonuclear electron neutrinos in Model 2.1 (including the 
normal mode) is shown in Figure As one can see, the total flux is quite large and 
significant up to E Uc rs 1 MeV. The flux from the massive component is cut off at E Ve rs 
0.1 MeV, which is far below the neutrino energy threshold at gallium experiments such as 
SAGE and Gallex-GNO {E threshold w 0.233MeV (221)- The normal mode thermonuclear flux 
is smooth and broad due to the larger SFR at low redshift. However, the flux is overshadowed 
by several orders of magnitude by solar neutrinos from the pp chains and the CNO cycle. 
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Figure 7: Fluxes from the massive modes of Model 2.2a (dashed red), 2.2ae (dashed brown), 
and 2.2b (dotted black). 



Reaction 


Endpoint 


% 


p+p^ e+ + v e 


0.42 


10 


13 N ^ 13 C + e + + „ e 


1.20 


45 


15Q _> is N + e + + ^ 


1.73 


45 



Table 3: These are the thermonuclear reactions considered here. The endpoint is indicated 
in MeV, and the percentage refers to the amount of the total thermonuclear neutrino flux 
due to each reaction mechanism. 



Since experiments sensitive to this energy range, like Borexino and LENS, do not have the 
capability to resolve any directional information about incoming neutrinos, the only current 
possibility for distinguishing SRN's from solar neutrinos is by spectral shape. 

5 Detection 

The detection of SRN's is inhibited mainly by difficulties excluding background events which 
include solar neutrinos, atmospheric neutrinos and antineutrinos, and antineutrinos from 
nuclear reactors. Cosmic ray muons also produce events that mimic the expected signal. 

The solar neutrino flux at Earth is larger than the expected flux from SRN's by several 
orders of magnitude for E v < 19 MeV [TT]. However, neutrinos, rather than antineutrinos, 
are produced in the thermonuclear reactions in the sun and have a smaller cross section for 
detection by about two orders of magnitude. This, with the directional information from 
recoil electrons in the detector, allows this background to be excluded at SK, KamLAND, 
and SNO. But decays of spalled nuclei from cosmic ray muons constitute an unavoidable 
background in this range. 

Atmospheric electron antineutrinos are cause for concern above ~ 8 MeV, but the flux 
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Figure 8: Total fluxes for Model 2.1(solid blue), 2.2a (dashed red), and 2.2b (dotted black). 




Figure 9: Total fluxes for Model 2.1e (solid cyan) and 2.2ae (dashed brown). 



becomes larger than the expected flux of SRN's only for E v > 35 MeV J2|- An atmospheric 
muon neutrino can interact with a nucleus to form a muon, which will be invisible in Cerenkov 
detectors if its kinetic energy is below the Cerenkov radiation threshold of 53 MeV. This 
background is significant for 19 MeV< E u < 35 MeV, but it can be described by the Michel 
spectrum and was subtracted off to obtain the current upper limit for the flux of SRN's at 
SK of 1.2 cm-V 1 for E v > 19.3 MeV [23]. 

It has been pointed out recently [23] that if the background analysis from SK is coupled 
with the sensitivity to electron neutrinos at SNO it will be possible to reduce the upper 
limit on the flux of electron neutrinos. SNO should be sensitive to a flux of 6 cm _2 s~ 1 in 
the range 22.5 MeV< E Ve < 32.5 MeV, which is an improvement on the Mont Blanc limit by 
three orders of magnitude. 
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Figure 10: Total flux of thermonuclear neutrinos from Model 2.1. 
In Figure ^2 we show the observable flux 

r°° dF 

F(E thresh )= / — dE (22) 

u ^thresh 

as a function of detector threshold energy. While the fluxes are quite appreciable at low 
threshold energies, they in fact remain relatively high at larger energies due to the large 
SFR associated with Model 0. Indeed, in all of our stage 2 models, our predicted flux above 
19.3 MeV already exceeds the current bound of 1.2 cm _1 s _1 from SuperK [231- The detailed 
flux predictions are given in Table 01 where we show the detectable flux for the viable energy 
windows at SK and SNO. Although SRN's will likely not be seen at SNO given these flux 
levels, in many of our models the SK bound is saturated by the expected flux, indicating that 
SRN's may be observed in the near future. This is in agreement with previous arguments 
made in based on simplified evolution models as well as arguments based on SN1987A 

Despite the large fluxes displayed in Table 0] relative to the SK limit [231 > one can n °t 
conclude that the stage 2 models considered have already been excluded by experiment. 
There are of course many uncertainties built into our chemical evolution models as well as 
uncertainties in the adopted neutrino physics. For example, one of the differences between 
our stage 1 and stage 2 models, is our choice of the IMF. In stage 1 models, the slope 
of the IMF was fixed at 1.7 whereas in stage 2 models, it is fixed at 1.3. The impact of 
this difference lies in the strong suppression of massive stars, i.e. the precursors of neutrino 
producing supernovae. As one can see from the table, the steeper IMF sufficiently suppresses 
the neutrino flux to satisfy the SK bound. 

One can also see from Table EJthat the dominant contribution to the integrated flux above 
the SK threshold comes from normal mode stars. Even the more extreme model 2.1e only 
contributes 0.1 cm -2 s" 1 to the flux above 19.3 MeV. However, these fluxes are very sensitive 
to our assumed average neutrino energy. Recall our adopted value for E Pe is 15.3 MeV. In 
Fig. H2J we show the sensitivity of the flux above 19.3 MeV (F(19.3)) to the average neutrino 



16 




Figure 11: Detectable fluxes from Model 2.1 with (dashed blue) and without (solid blue) 
oscillations as a function of neutrino energy threshold. 

energy. In order to satisfy the SK limit of 1.2 cm~ 2 s~ 1 , we would have to lower (E Pe ) to 13.3 
MeV. This is fully consistent with the range of neutrino energies in supernova models |19j . 
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Figure 12: The integrated flux above 19.3 MeV in Model 2.1 as a function of the average 
neutrino energy. 

Ill-understood backgrounds and large detector energy thresholds make detection of SRN's 
difficult. Because the SRN flux peaks at or below 1 MeV, there is little hope that the peak 
will be probed with existing experiments. However, being able to exclude backgrounds at 
lower energies would greatly increase the probability for observation at SK or a similar 
experiment. 
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Model 


SK Flux 


SNO Flux 


1.0 


0.40 


0.11 


1.1 


0.40 


0.11 


1.2a 


0.40 


0.11 


1.2b 


0.40 


0.11 


2.0 


1.8 


0.47 


2.1 


1.8 


0.47 


2. lose 


3.2 


1.4 


2.1e 


1.9 


0.49 


2.2a 


1.9 


0.48 


2.2ae 


1.9 


0.49 


2.2b 


1.8 


0.47 


2.2bosc 


3.2 


1.4 



Table 4: Predicted fluxes in cm 2 s 1 in the models considered here. Results are given for 
electron antineutrinos with energies E Pe > 19.3MeV for SK and for electron neutrinos with 
22.5MeV< E Ve < 32.5MeV for SNO. 

6 Effects of Oscillations 

Due to our assumption of equipartition of energy among the neutrino species, oscillation 
will only affect the flux if different species emerge from the explosion with different aver- 
age energies. Indeed most supernova model calculations do show a hierarchy of neutrino 
energies, and because the energies of E Pfi and E Dt are generally higher than E Pe , the effect 
of oscillations will in general increase the observable flux. The effects of oscillations on the 
neutrino background were considered previously in [26 . 

Figure IT3~1 shows the total flux of electron antineutrinos from Model 2.1 both with and 
without oscillations, where we use the neutrino average energies discussed in section El and 
maximal mixing has been assumed. Although the total number of electron antineutrinos 
arriving at Earth is smaller due to oscillations by ~ 16%, the flux of neutrinos with energies 
greater than ~ 9.5 MeV is larger. The differential flux for Model 2.1 with and without 
oscillations is shown in Fig. ED Note that the solid curve here is identical to that in Fig. |HJ 
Here we clearly see the individual peaks due to the massive mode at < 1 MeV and normal 
mode at ~ 3 MeV. 

For E v > 19.3 MeV, oscillation effects would increase the observable flux at SK by as much 
as 78%. This effect is seen in Fig. ITT1 where one sees that the integrated flux with oscillations 
(dashed curve) exceeds the flux when oscillations are ignored for threshold energies greater 
than about 6 MeV. The effect of oscillations on the SK observable flux is seen in Table 0] 
for the models labeled 2. lose and 2.2bosc. Reconciling these fluxes with the SK limit would 
require a further drop in the average neutrino energy or a tightening of the assumed neutrino 
energy hierarchy. Similarly, the flux of electron neutrinos potentially observable at SNO is 
increased by almost 200% in the energy window 22.5MeV< E < 32.5MeV. Although our flux 
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calculated with oscillations is still less than 1/4 that necessary to approach the projected 
SNO sensitivity [21] , with a better understanding of backgrounds the prospects for detection 
in the near future are encouraging. 




5 10 15 20 

E v (MeV) 



Figure 13: Total fluxes of anti-electron neutrios from Model 2.1 with (dashed blue) and 
without (solid blue) oscillations. 



7 Conclusions 

We have considered several scenarios for star formation which reproduce the observed chem- 
ical abundances and SFR for z < 6 and reionize the universe at high redshift. Each model of 
star formation here consists of a normal mode coupled to a Population III mode of massive 
star formation at high redshift. We examined the SRN flux from the core collapse supernova 
explosion as well as from the thermonuclear burning stage. 

Because the massive mode of star formation is so brief and takes place at high redshift, 
the corresponding electron anti-neutrino fluxes peak at E v < 1 MeV. Thus despite the large 
fluxes produced by the massive mode, these low energy neutrinos will be difficult to detect. In 
contrast, the normal mode of star formation, which dominates the flux at observable energies, 
is peaked at a somewhat higher energy and has a broad spectrum due to the production of 
stars at lower redshift. The neutrinos produced during thermonuclear burning are emitted 
with much lower energies, on average, but also exhibit the behavior of a sharp peak due to the 
massive mode added to a broader spectrum from the normal mode. Thermonuclear neutrinos 
are unlikely to be observed at experiments which cannot at least partially resolve the direction 
of the incoming particle because they are produced in the electron flavor state and the 
spectrum lies several orders of magnitude below that of solar neutrinos. Our calculated 
fluxes of SRN's from core collapse, however, saturate the SK bound of 1.2cm _2 s _1 for E v > 
19.3 MeV in all stage 2 models. Although there are uncertainties in the neutrino physics, 
such as the average energies at emission, the prospects for observation in the near future 
are good. We also examined the effect of oscillations by calculating the flux with maximal 



19 



mixing. With the accepted neutrino average energy hierarchy, (E Ue ) < (E„ e ) < (E Ux ), any 
oscillation will harden the high energy tail of both the v e and the v e spectra. 

Further refinement of the neutrino physics and measurement of the SFR out to higher 
redshift would allow for a more definite flux prediction. With decreased detector thresholds 
and increased background rejection, observation of the SRN flux will soon be possible. 
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